C 


NACA TN 3224 91G9 


| 


-~ — - 


p 
NATIONAL ADVISORY COMMITTEE 2333 
FOR AERONAUTICS B 


TECHNICAL NOTE 3224 


THEORETICAL INVESTIGATION OF THE EFFECTS UPON 
LIFT OF A GAP BETWEEN WING AN D BODY OF A 
SLENDER WING-BODY COMBINATION 
By Duane W. Dugan and Katsumi Hikido 


Ames Aeronautical Laporatory 
Moffett Field, Calif. 


WASHINGTON 
August 1954 
Ai ELE 
20 WAP 
> AY pa da CR ree i : " 
Ceu uc lonas 
RON 


TECH LIBRARY KAFB, NM 


NATIONAL ADVISORY COMMITTEE FOR — ANE 


O44 


TECHNICAL NOTE 3224 


THEORETICAL INVESTIGATION OF THE EFFECTS UPON 
LIFT OF A GAP BETWEEN WING AND BODY OF A 
SLENDER WING-BODY COMBINATION 


By Duane W. Dugan and Katsumi Hikido 
SUMMARY 


Slender-body theory has been applied to the study of the effects 
upon lift produced by the presence of a gap between wing and cylindrical 
body of a slender wing-body combination. Two conditions were studied, 
one in which both wing and body had the same angle of attack, and the 
other in which only the deflected wing had an angle of incidence to the 
free stream. The lift for the case of combined angle of attack and wing 
deflection can be found by superposition. 


The theory predicts large losses in lift even for minute gap widths; 
it is anticipated, however, that the effects of viscosity and possibly 
of compressibility not considered in the theory will serve to reduce such 
losses in practice in the case of very small gap widths. The loss in 
lift effectiveness due to gap effects is more severe when both wing and 
body are at an angle of attack than when only the wing has incidence 
relative to the free stream. For the wing-body combination exclusive 
of the nose and afterbody, the gap effects, expressed in terms of per- 
cent loss in lift, are more pronounced for larger ratios of body radius 
to wing semispan; and the ratio of the lift obtained from wing deflection 
to that due to angle of attack increases with increasing gap width. The 
effect of the lift of the nose is to reduce the percent loss in lift due 
to gap in the angle-of-attack case from that predicted for the portion 
of the configuration considered above. As a consequence, the effective- 
ness of the wing as a control surface may increase or decrease with 
increasing gap width accordingly as the ratio of body radius to wing 
semispan is small or large. 


INTRODUCTION 


In connection with the use of all-movable lifting surfaces in 
missile design, there arises the practical problem of the effects upon 
the aerodynamic characteristics of the missile caused by the presence of 
a gap between the wing panels and the fuselage. When the fuselage is 
cylindrical, the gap is unavoidable for two reasons: In the first place, 
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a clearance between the movable wing panels and the body is required 
from mechanical considerations, and secondly, a space between the wing 
panels and the curved surface of the fuselage is created by the deflec- 
tion of the wing with respect to the body. This latter gap distance 
varies in the chordwise direction. For the usual ratios of wing chord 
to body radius encountered in missile design, and for small angles of 
wing deflection, this variable portion of the gap width will be very 
small compared to the constant gap present at zero deflection. In the 
present theoretical treatment of the problem, only small angles of wing 
deflection are considered, and therefore the gap is assumed to have 
constant width. 


The purpose of this report is to evaluate the effects of gap upon 
the lift of a typical slender wing-body combination having an all-movable 
wing. Two basic problems are considered: Problem one deals with the 
lifting characteristics of a wing-body combination inclined at a small 
angle of attack with respect to the free stream, the wing panels having 
zero deflection with respect to the body; problem two is concerned with 
the same configuration, but the body is considered to be at zero angle 
of attack and the wing is deflected to some small angle with reference 
to the body axis. By the principle of superposition, the lift of the 
wing-body combination having both angle of attack and wing deflection can 
be found. l 


Although the first of the problems described above could be treated 
by the method given in reference 1, where the effects of a gap in the 
middle part of a wing at subsonic speeds are investigated, the second 
problem appears less amenable to solution by such methods. More useful 
for the present purpose, it is believed, is the procedure employed in 
reference 2, in which slender-body theory is applied to the analysis of 
the aerodynamic characteristics of a family of wing-body combinations. 
Accordingly, both problems are treated in much the same manner as that 
given in the latter reference. 


While the present report was in preparation, the results of other 
investigations of gap effects were published. Bleviss and Struble (ref. 
3) suggested a method of estimating effects of streamwise gaps upon the 
lift of wing-body combinations which included replacing the body by an 
infinite wall. Mirels (ref. 4) obtained the slender-body result for 
problem one of the present analysis, but did not obtain in closed form 
the solution of the second problem. 


In common with the above-mentioned investigations, the present 
analysis makes use of several simplifying assumptions. For example, the 
wing is considered to be vanishingly thin, and the fluid medium is 
assumed to have zero viscosity and to be free from shock waves. Con- 
sequently, the usual side-edge condition of infinite velocity and zero 
lift is applied at the edge of each wing panel adjacent to the gap. It 
is also assumed that the lift is zero across the gap. Only small angles 
of attack and wing deflection are permitted, and the plane of the wing 
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is considered to include the center line of the body. Finally, the wing- 
gap-body combination is regarded as being slender in the sense that 
lateral dimensions are small in comparison to longitudinal (free-stream) 
dimensions. 


For gap widths of the order of a few boundary-layer thicknesses, the 
neglect of viscosity and compressibility will certainly produce unrealis- 
tic results. For example, subsonic theory (ref. 1) predicts large losses 
in lift for even infinitesimal gaps in the middle part of a wing. It is 
interesting to note, however, that in an early experiment treating the 
configuration of reference 1, Munk and Cario (ref. 5) obtained results 
which indicate appreciable decreases in lift-curve slope at zero angle 
of attack for gap widths which, though not "infinitesimal," were small 
in terms of the wing chord. It remains for experimental investigation 
at supersonic speeds to define the range of gap dimensions and of angles 
of attack for which the present inviscid theory is reasonably valid. 


GENERAL ANALYSIS 


If it is assumed that for slender bodies (including wings) the 
velocity gradient in the free-stream direction ?y, is small compared 
with the velocity gradients in the y and z directions, then the well- 
known Prandtl linearized differential equation for the perturbation 
velocity potential € of a compressible flow in three dimensions 


is closely approximated by the more simple twa-dimensional form 
Pyy + Paz = 0 


which permits the use of conformal mapping. 
The same simplification can be made if Mo, 
the free-stream Mach number, is unity, or 
very close to unity, and the velocity gra- 
dient in the free-stream direction is not 
excessively large. 


The wing-body combination chosen for 
consideration is shown in sketch (a). 
Symbols are defined in Appendix A. The 
body consists of a cylinder of circular 
eross section of constant radius and an 
arbitrarily pointed nose. The wing con- 
sists of two identical right triangular 


Sketch (a) 
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flat plates? of vanishing thickness placed adjacent to and on opposite 
sides of the constant cross-section portion of the body. The distance 
from the body axis to the inner edge of.each wing panel is designated 

t; the gap is defined as the distance t-ro, where ro is the radius of 
the cylinder. The distance from the body axis to the outermost point on 
the wing panel is called so; the symbol s refers to the distance from 
the body axis to a point on the leading edge of each wing panel. The 
slope of the leading edge may be designated as ds/àx, or m = tan €, 
where «€ is the vertex angle of the wing panels. The origin of the 
coordinate system is taken at the point of intersection of the body axis 
and the line joining the apexes of the wing panels. 


We introduce the complex variable & = y+ iz = peł, The cross 
section of the wing and body in the E plane is shown in the upper part 
of sketch (b). The Joukowsky transformation 


2 
r 10 
E = E pS E, = Yi + iz, = pje * (2) 


E 


maps the given cross section onto the real axis of the §, plane, as 
shown in the lower part of sketch (b). 


iz 
E plane 
ro 
-8 , =t t S x 
-B31 “ti arı I4 tı B1 ^. 
Sketch (b) 


lThe results of this analysis are, however, applicable to con- 
figurations in which the leading edge is not a straight line, provided, 
of course, that the assumptions of slender body theory &re not violated. 
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Then, for 2, =0, 


ro 
maS a an 
yı = 2rọo cos 0 = Ly yi^ € n? 
B 
81 = S + —- (3) 
B 
2 
ro 
ti = t + — 
t t 
rı = 2rg 
and 
NC NGA 
y = Xi Ji “1 y; r 
2 
2 2 
Yu - Vy -ri^- n 
y = = nír M 
2 
yeu y? e r? 


The induced velocities in the two planes are related by 


v - iw = (vi-iwi) ae (5) 
as 
from which 
: ro Y r = 
v= v |i (2) cos 20 | + (2) sin 26 
p p 
ro s Io a . 
w= "ww |l- 7 cos 20 | - v4 7 sin 20 (6) 
| Eoss Yo \2 
Vp =( Vy cos O + Wy sin 8) i= >) | * (2) sin 0 


Also, Laplace's equation must be satisfied in the E, plane, so that 


Pyaya + Pzyz = O (7) 
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The general solution to equation (7) that gives the vertical induced 
velocity w4 on the yı axis due to the jump in the value of the lateral 
velocity Avy, across the yı axis can be written 


l 
wan) 2- Bf a (8) 


ex -81 Ji 7 Ya 


In the present problems, w,(y,) is known from boundary conditions, but 
Av; is unknown (except in the gap). Equation (8) is thus an integral 
equation which will be solved to find Av, in each of the two problems. 


PROBLEM 1 - WING-BODY COMBINATION AT SMALL ANGLE 
OF ATTACK, WING AT ZERO ANGLE OF DEFLECTION 


Boundary Conditions 


As was mentioned in the introduction, it is assumed that no lift is 
carried in the gap; hence the circulation and vorticity in the gap are 
taken to be zero. The boundary conditions which obtain in this problem 
are therefore 


(i) w= -Voa 02 rg, 0508 2x 

(ii) Av = 0 zZ=0,ro <y'< + 
(111) w= -Voa z =0, £f < y< gl (9) 
(iv) v=0, w=0 p=w, OS 6 < 2x 

(v) AP=0 z = 0, y 2 8 


In the 5, plane these boundary conditions become 


Vo 
(1) eS Zi = 0, 0€ y? € r$? 
(ii) Avy, =0 23 0, Ya^ < nE < tf 
Vom yı 
(iii) w = - — [1 21 = 0, t4* « yi? < 81° (10) 
2 yi? =p, 2 
(iv) va =0, w= 0 0,29, 0€ 04 < On 
(v) AQ, =0 zy = 0, y? > s,f 
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Solution of the Integral Equation 


The general solution expressed by equation (8) can be written 


1 Av 1 l Av 
wily1) = - — IL ay, - — —_£- ay, - 
: ex “Sz M MA fo 2n -t1 yi = yo 


— ——— dy. M m —À dys ie AN dyz 
et Jary Yi ~ Ye et dp, Yi 7 Ye ex Jt. Yi - Ye 


where 


AViy = AV between ti and +84 
Aig = Avı between trj, and ttj 


AV4g = Av, between -rı and rı 


Since Avig = O, the second and fourth integrals of equation (11) vanish. 
Also, since the wing-body combination is laterally symmetric, 


Avily1) = -Avi(-y1) (12) 


and equation (11) can be written 


mtm) = - > [ ¡PE AE [7 AW an (13) 
23 Jg "ni 7 "5 2x t2 Ny 7 No 
where 
4 = yi? 
Me = Ya” 


Inversion of equation (13) will be accomplished by the application 
of the procedure presented in reference 2. First consider the region 


Voa 
os n, < r?. Here wa =- ra and equation (13) may be written 


2 2 

Voa m: Pan Av si (02) T ES [ Avi (n2) à 
V c. e Pa 
bi Ta 5 No en O 


cw A AA — — A A AA —— a on 


i ——— — —M— € — ——— ——— ee 
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Equation (14) is now of the general form 


b Avi (n4) 
f(n, ) = ej 1, - Ty dn, (15) 


of which the inversion is 


Avı (A) = 2 Fi ^r; (n,) dn, 2 £n) / (5-0, ) (n, -8) à 
tf (b-A) (A-a) a a A =a 


(16) 


which for Avi(a) = O reduces to 


i f b- 
ann) = 2 fA fU uia [78 (27) 
4 


In the present problem, Avı (0O) = O from symmetry, and the inversion 
of equation (14) is 


81% Ay a a) 
Apl) = - Voa wh — |. Es ber VÉ 
2, y x fr? m t? os 


(18) 


This expression for Avy, will be used subsequently in solving for the 
lifting properties of the body. To find Ay it is necessary first to 
substitute equation (18) into equation (13) and reverse the order of 


V LI 
integration. Using wi(n, ) ae l + e for t? < n, < 51° 
T. 7r1 
1 


we obtain, after simplification, 


2 2 
81% Avalon Y /nm.-T1 
1 1 


Inversion of equation (19) will follow the procedure given by 
equations (15) and (16). Thus 


2N 
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2/7 
EAN e Ce 
AViy = 


Ta (12 7n, ) (n, -t15 ) (n, 11) 
- nai Y (gt) a 


t1? Nh - No 


(20) 


or 


2 
Ava = DO A { - Voa ^ len, - (5*9) ) 
2,/(812-9,)(n,=12)(n, 112) 


(21) 


where 


1 


Equation (21) is not unique since it contains an integral involving 
Aviy. However, the integral is merely a constant, which can be evaluated 
from the assumption that A®, is zero at yy" = ti. Knowing Avı, we 
can express A, as 


Ji 
AQ, = [ Avi dyy (22) 
5 al 
or, in terms of i a 
"a Av 
AT, = [ E an (23) 
2 1 
81 2 "a 
Then 
t^ ^v 
Am) =0= | H an, (2h) 
81 D Ma, 


Substituting equation (21) into equation (2h) yields 


ese: e me _— a, — ee n a A eG A M aa I, AI en AA ER RR a A A Án — 
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2 
A A ucc NE osi 
2 2 2 l 
B S3 
* X (8:2-n,) (n47t:2) (04 23?) JV (817 =n, ) (n, =t13) (n, 7037) 
(25) 
Performing the indicated integrations? and solving for G gives 
G = 1Voa | (s17-r3*) Ee l (5,2+t,2-21,2) (26) 
K(k) 2 
where 
81^ -r3* 


An expression for Avıy can now be found by substituting equation (26) 
into equation (21). The result is 


AT = 2 Ani Vou (mors?) EOD. ~ (9, =r) | (27) 


Jf (837 =n, ) (n, -t3?) (n, -115) E) 


Span Loading on Wing Panel 


The span loading is proportional to the total circulation 
T = AP - AP. From equations (23) and (27) we obtain 


APiy = 2V ga. l / 941? -rQ* z(y 1) k) (28) 


where 
Z(Y K) = E(¥,,k) = — F (y, ,k) 
and 


Y. = sin”? (a 
81% -t1 
(The function Z is tabulated in ref. 7.) 


“The elliptic integral of the third kind encountered here, as well 
as those which occur subsequently, were evaluated by use of reference 6. 
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If t4 is set equal to r, (condition for a slender wing on an 
infinite cylinder with no gap) equation (33) becomes 


(A93). or, = 2Vya,/ 837 -y3* 


which when transformed by equations (3) may be written in the yz plane 


Tot rot | 
(AP) par. = Poo | e? € + — -yP (1+ - (29) 


Equation (29) agrees with the expression for the velocity potential for 
a slender wing on a semi-infinite cylinder given in reference 8. 


At the leading edge of the wing, yı = 81, and (AP 1. = 0 as it 


should. Therefore, Ty can be given by equation (28) with 8/7 substituted 
for s. However, it is more convenient to express the value of the 
circulation in terms of the value of the semispan of the wing for zero 

gap Sg*. Then, since Sọ = So* + g, where gœ is the gap distance, the 
expression for span loading is 


a 2 | 
an \ B r 
= = Volgy = soa e + £) |: - E l Z(¥, ,k) (30) 


where Y, and k are evaluated for 58] = Sọ- 
Pressure Distribution Over Wing Panel 
The loading coefficient Apta can be found from the equation 


Ap _p¿M_2 JAP 2 AP ds, ds 
a "Vo Vo Ox Vg Os, ds dx 


First, the differentiation of A?,,, with respect to si is carried out 
and gives, after simplification, 


JAP 7 81 y 81^ -ri* E(k) (y1? -t1*) (y1*-r1?) 
——= = 240 1 - —— |[Z(Ouk) + [| 
981 81? - tf K(k) (837? -y32) (s3* -r3*) 


(31) 


— — o A A A cl dn ne c nc md -— —— M M —— MP A A —— pn  —M P —M M 2 a AA Pa qu TT A A A eS MM — -n e ——— 


==> - — — ee — — mL aA LL. 
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Since ds,/ds = 1 - roọ?/s2, and for triangular wings ds/dx =m, there 
results for the loading coefficient for the wing 


oy” tn O) E) tem 


e -12 y) GE aot) 
yea (PoP Pet) ide 


In the case of zero gap (t = ro), equation (32) reduces to the 
result given in reference 8 for the loading coefficient for a triangular 
wing on & semi-infinite cylindrical body. 


Lift of Wing Panels 


The lift of the wing panels may be found by integrating the 
expression for the spsn loading; that is, 


OR- Lo mra t (33) 


The integration is more easily carried out in terms of yı. Equation 
(33) then becomes 


i 


(3 ra d BEC C + JE) dyi (34) 
1 -rif 


Bübstituting in equation (34) the expression for > CA op obtained 
from equation (28) with sı = soi gives 


€ : = ri $ (801 217) k- 2 560) I [2+ Ao (5/3) - qu? + 
2t, 3012-1? (2082) - K(k1) + S E - ste) Cs) | (35) 
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where 
Ao Wri) = 24] 80a) - x03) | Piia) + 803 )) 
Y = sin? [STi y 2Iig ri Pt 
s42 ta ta [| g42-r? 
ki! = JA - ky? 


and the moduli and arguments are evaluated at the wing trailing edge 
where 81 = So. (The function Ag is tabulated in refs. 6 and 9.) 


Equation (35) gives the lift of two pointed low-aspect-ratio wings 
in the presence of a slender body which is cylindrical in the region 
between the apexes and trailing edges of the wing panels and lies 
between and at a distance t - ro from the panels, all at an angle of 
attack. 


If ti = rj = 0, equation (35) reduces to (L/qa)y = 218,2 or, for 
the two triangular wing panels now joined (triangular wing) 


218 


Cy, g Aa ( 36) 


m a 


which is the well-known result for the lift of & slender wing. 
Span Loading on Body 


Equation (18) expresses the vorticity on the body in terms of the 
vorticity on the wing. Substituting equation (27) into equation (18) 
gives 


817 -r4*)d 
Avr p (n,) = = Vou Ny la I i A Mrs 000 - 
ri?-n t^ (nima) / (81? -n,) (n,-t:2) 
a so (^r — ——— MAE. 
ieri K(x) 


(nina) / (812-14) (n,-t17) 


A A A AA KA AAA ee un 
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which -becomes 


2) EQ 


(ri*-71,) + (81-1 
Avipln1) = Woa / = - o 5:2; (38) 
a JA 83? -n,) (37 -n,) 


m |e 


Since A, = O in the gap, the jump in potential across the real axis 
between -rı and rj can be given by 


y Av 
Aage | 2 dn, (39) 


From equations (38) and (39), it follows that 


1 
2 V ya F s1?-yi? sin Yo - ya ri? Z(Yo,k) - > J neve | 


. 2 2 
Y = sin”+ I1 731 Y 
= 4,2 yn 2 
a "JL 


ri, equation (lO) reduces to 


we B um ub 2 v2 
(Ag). or = eVoa RE yi 5 J ri^ yz | 


or, in physical terms, to 


2\2 
(APB) Ln = 2Vya | /s* € + 29) - hy? - /r,?-y? 


(o) 


where 


When ty 
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The span loading on the body due to the presence of the lifting wing 
can now be determined. First we shall consider the loading induced on 
the body aft of the apex of the wing 
(see sketch (c)). The value of 
(A915) 7 is found by setting s, = tx 


in equation (40), which becomes 


(A9, rx = Vya R rj? -ya 


For APig at the trailing edge, the 
value of So, is used in place of 8; 
in equation (40). Therefore, the span 
loading on the body is 


Sketch (c) 
E = haa E Soi--yi^ sin Ya - ,jSoi?-ri* Z(W45,k) - y zm | (41) 


or since the value of the span loading on the nose in Nga ,/ ro -yF 


according to Munk's airship theory, the total span loading on the body 
is, in terms of the geometric properties, 


dL : 
“a 2\2 2 2 
dy E d To y To 
E o” | | Bo Bo- E Bo* ; eus 
B 


NC 2 
Tu : 
Bo Bo 
where Y, and k are evaluated using 81 = 801. 


From equations (28), (41), and (42), the span loading for the wing- 
body combination is plotted for several ratios of gap to semispan in 
figure l(a). No lift is induced on the body aft of the trailing edge of 
the wing, according to slender-body theory. 


“mn NENNEN "i |J ann NT mó. a e 
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Pressure Distribution on Body 


The expression for the loading coefficient Ap/q on the body is 
obtained in the same manner as for the wing in a preceding section. The 
result is 


DEOL 


(43) 


If t = rg, equation (43) reduces to the result given in reference 8 for 
the lifting pressure coefficient on an infinite cylindrical body in the 
presence of a lifting low-aspect-ratio wing. 


The distribution of lifting pressures on & spanwise strip for the 
present wing-body combination is shown for several ratios of gap to 
semispan in figure 2(a). For comparison, the pressure distribution of 
the wing alone is also shown in the same figure. 


Lift of Body 


The lift induced on the body by the wing can be found by the 
formula 


To at rı 
wf Raf B an (a) 
-rọ Y o Gy 


with dIg/dyi given in equation (41) above. The final result can be 


written 


(=), = D + (5o. rs) |a Eus = SANE S Ao a). - uj? - 


2ta sor? | Ba) - K(k1) + us E - 2 SH fca) ) 


(45) 


where Y, k, and kj are evaluated for sı = SoL- 


3N 
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For the lift of the wing-body combination (exclusive of the nose), 
equation (45) is added to equation (35). This yields the simple equation 


€3 à ex [ (ons) + E - 2 m] (46) 
WB 


When written in terms of the physical dimensions, equation (46) agrees 
with equation (2) of reference }. 


The total lift of the body and of the wing-body combination can be 
obtained by simply adding the lift of the nose (which has the value 


x 
z r? = 2nrof) to the respective expressions for lift. 


For the case in which the gap is zero (t, = ri), equation (46) 
reduces to the result given for this case in reference 8. 


Equations (35), (45), and (46) have been evaluated for several gap 
semispan ratios; the results, given as ratios of lift in the presence of 
a gap to the lift when the gap is zero, are shown in figures 3(a) and 


la). 


PROBLEM 2 - BODY AT ZERO ANGLE OF ATTACK TO FREE STREAM, 
WING DEFLECTED AT SMALL ANGLE WITH RESPECT TO BODY 


The second problem is complicated by the fact that the total gap 
between the wing panel and the cylindrical body is the sum of two gaps, 
one of which is the constant clearance required for mechanical reasons, 
and the other is the gap caused by wing deflection. The latter gap 
distance varies from apex to hinge line, and from hinge line to trailing 
edge of each wing panel, and also changes in magnitude with the angle of 
deflection. However, for the small deflections considered here (say of 
the order of 5°), the gap created by deflecting the wing will be small 
for the usual ratios of body diameter to wing chord encountered in prac- 
tice, and the total gap distance will be taken as a distance independent 
of wing deflection angle. 


Boundary Conditions 


Certain boundary conditions for this problem differ from those of 
the previous problem, although the basic conditions, such as the zero 


ae een aie ee em a MP y mos a 8 m mme S nh s A t RI a gio t pana mee ae n Ne ee a a e m m m. a 
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value for the jump in potential at the leading edge of the wing panel, 
apply here as well. Stated first in the physical yz plane (see sketch 
(b)), the salient boundary conditions for this problem are 


(i) vr =0 p=xro, OL 0% 2x 
(ii) Av =0 z = 0; Boo «y^ e 
(iii) w= -Vod z=0, tf <y <? (47) 


(iv) v=0,w=0 p=0, 0[9<?2x 
y 2s 


Transforming equations (47) to the £, plane by means of equations (6) 
gives 


li 
O 


(v) AQ 


(i) w =0 zı 2-0, 0« y? <r 

(ii) Avy =0 Zi = 0, ry? < yj? < t4? 
(iii) w =~ Lm R + === tif < yy? < s; (48) 
(iv) vy = 0, wy = 0 Duo P =m, OL 9 Sea 

(v) AP, = 0 yi* 281 


Solution of the Integral Equation 


Commencing with equation (13) we shall obtain a solution for Avi, 
the jump in lateral velocity along the real axis. In the region 
O< n.-< rı? the vertical component of velocity is zero so that equation 
1 
(13) can be written 


2 
81 Ay ri Ay 
x —— an = - — ài; 0< n, € ri? (49) 
x tf 74 7 No er O T4 7 No 


Applying the inversion formulas given by equations (15) and (17) gives 


e 
[n 71 d 
AT pl M2) a E E " > 
n? JT“ Jo ne-m 


2 
1-0, 1 Aval na) 
—— cmd Ut 
"na t2 M7 Ma 


(50) 


- e e vex ae) eee e ee m 
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or after integration with respect to US 


( ) tub No 81 Avy (n3) n47r3* 
ratte) “= fata ae ng OS 
2 Vt? 2 3 3 


Next, substitution of equation (51) into equation (13) will eliminate 
Avipg and provide an integral equation for the vorticity across the wing. 


Thus 


2 + 2 
y ( ae a tele rı . dn. No B1 Ava (Na) a-r" à 
a(n) = - 5 cep re peu eg 
O Ny i No Ii “Neo tn No - Ne Ns 


2 
83 AV, 
ti^ Ny No 


Integrating the double integral first with respect to n, We find that 
for ti^ < n, < 817, equation (52) reduces to 


2 
2. 81 2 
Nr 1 Armin.) Ma-ri 
1 ña 2x 2 My” Ma Ns » 


Equation (53) can be inverted by the formulas given by equations (15) 
and. (16). Thus is obtained 


A 


Aviy = A | c + 


x, (812-1, ) (n, -t12) (n, -118) 


(8317-15) (n5 -t37) (n5 2737) 
Np 


j w1 (95) 


dn (5h) 
tf "a Ei Yo = 


ta ee re A a PP XX e aa, 
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where 


81 
t4 


As in problem 1, the value of the constant C can be determined from 
the condition that A,(t1) = 0. Since . 


and from boundary conditions 


Vod Mo 
jr fix J 
2 2 ( 35) 


it follows from equation (54) that 


2 
tí dn, 


t an Vod 
54> (8,20. ) (n, -t32) (n, 042) 2 Jg? 2 2 2 
8i “9/1 7*1 AT, “Pa (81 -5, )(n, -t1 )(n, -r1 ) 


(81° -n3 ne -t Mn, -11?) 
——M———— 4. 


Mo 


t^ dm, 
ri — f(81*-n,) (nz -t37) 
8 e Ty “Ns 


(55) 
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from which 
Vo /s1?-ri? t dn, 
y K(k J 2 
C9 Men Ce-a, (n, tx?) (n, 23?) 


81-92) (no-ta*) (n5 -r17 ) 


X 
: -5 (en, - 81^ - t°) 
e 


(56) 


This expression is left in this form since further integration is not 
advantageous at this time. 


Span Loading on Wing Panel 


The circulation A9,(y1) will be obtained by integrating the 
expression for Avyy given by equation (54); that is 


2 : 
Ji Av 
AP (yi) = Fi = dn; Mm = a” (57) 


81% 2 fh 


By a suitable choice of the order of integration, equation (57) can be 
evaluated to give 


ri? 


APh ya) = Voð {2 K(k, ) | Z(V. ,k) = nitro ska) | + 


1 


Jaén? ZYK) | L+ ned |} (58) 


(51 -y,7)t,* 


(512-415) y1* 


where 


" — om. -— — = — 
ee me apa com an LIP P A A — a ca - A r 
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In the case of zero gap (t, = ri), equation (58) degenerates into 


Vod r 
A34 mr, ) = = x 91? -yf E + = cos”+ = | + 


T Si 
r 
/ sy -rj? + x 81? -yi* 
yiln - - 
2 2 Ii 2 
SI XE. cmo 81 “Y1 
yi 
2 2 2 2 
JJ 81 -rı + 81 -y1 
1 GA —_—_ (59) 
2 2 2 
af 81 713 - S1 -y1 


Equation (59) agrees with an unpublished result obtained by Gaynor J. 
Adams, who solved a similar problem of a slender wing-body combination 
having the wing deflected - no gap being considered. 


Since AP= O at the leading edge, the expression for the span 
loading can be written 


2 
aL To 
d 8 
E nc que. TEL NL K(k,) + 
hagso x So So tf ro 
Bo" Bom 
2 
E 
£m fa. + aoti] Z(Y, 5k) - 
Bo 
2 2 
s T 
(+ E )x(%) Z(Va,k1) (60) 
Y NES Be? 
O O 


where the moduli and arguments of the elliptic integrals are evaluated 
for 81 = 801- 
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Pressure Distribution on Wing Panels 


Following the procedure outlined in problem 1, page ll, we solve 
for the loading coefficient of the wing panels. Differentiation of 
equation (58) with respect to sí yields 


SA, |: 


p - m ge t? ([ a it mon? E PAN w + 


er," 


= x(n)  - 2t; [K(%r) - Ba) 1 ) Ata + 


2-t12) (y1*-r19) 


(83° -y1 ) (833? -1,8) 


(yı (61) 


and the expression for the loading coefficient becomes 


AN _2m_ (s^-ro*)t? PPE) f(y trot RA 
qd /y x s(s*-t®)(s*t®-rot) | yt(s®-r,7) / (s?y? -rg*) (s? y) 
Org^t | t 
Z(Y,,k) ({2- T "ir wes AE 


8 


x(s*-ro^) E " Ao] } = aly emo | KO) = EG) |) 


(62) 


— —À 
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For the limiting case of zero gap, equation (62) degenerates to 


^p _ oa [77 G) | ES (63) 
Bea TEOTEO 


Lift of Wing Panels 


The solution for the lift of the wing panels in this problem follows 
the same pattern as adopted in the first problem. In the €&, plane 


L 2 Bol yi 
LN -£ (A9) (4 + ==) dy (64) 
( Hu Vo J. | TE / yi? -r12 L 


Substitute the expression for Amp obtained from equation (58) into 
equation (64) and integrate and there results 


2 
E(k) 2r1?K( ky) 
€ 3h aG e IE i i T d "n ui v T c fea 
ety / 801? -r3? | BCs) - m 1 + Ao(U ,14 ) 2 rj^K(k;) P 


E tty /891?-rj? 


2 
2 - E(k, )K(ki) - : rj? (65) 


where the moduli and arguments of the elliptic integrals are evaluated 
using 81 = 801° 
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Span Loading on the Body 


Equation (51) provides an expression for the jump in lateral 
velocity Avig along the real axis in the transformed plane between 


O and ri^ in terms of Avy between tı? and 58,2. Combining equations 
(51) and (54) gives the expression 


2 
Avy p(n, ) = E hes de dna C + 
pa 8? (nina) /(8317 7n.) (n, -615) 


ps (ei ng) (netz (ei) , (66) 
t12 (na-na) = E 


An expression for APig can now be found by substituting equation (66) 
into 


Using the expression for C given in equation (56) and performing 
the integrations gives 


V 5 = 8 =~ t= - 
AM, m 2K( kz) l rit-Yi^ 1” Saya 81 Wr Ya” yi^ 
n 1 ea? Ly]? 81^ -yf 
ZU z(a,k) = alok) | - x ma? 
E 2 “YL 43 51 1 “JL + 
1 


x E + no(t) | | J 81 -ya^ sin Va - /fs*-ri* Zítka 1) | 


(67) 


= a ʻa -—- - ——-- —— an — —— — —À— —— À MÀ M M P À€— MÀ €  — M t€ m I —— e ae 
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where 


( 817 -rf ) yi* 


( 812 -yi* jr? 

With 8, = Sol; equation (67) gives the value of APg at the 
trailing edge of the wing. If the body extends aft of the trailing edge 
of the wing, the present boundary conditions cannot properly be applied 


to find the loading on the afterbody. No attempt to determine such load- 
ing will be made here. 


Since (APip) ig is zero, 


2 
rı 
ty Z(V o, K) = y1Z (Y, ,k1) - T pri -y1? T 


t[1- Ay(Y,k1)] F 8o1^-yi^ sin Ya - / 891? -rj* 20.3) | 
(68) 


where 801, the value of s, at the trailing edge, is used throughout. 
For the condition of zero gap, the expression for span loading can be 
shown to be 


al; 
=) P| a( [nos 2 y - fri -y3 $ + 2 /801% -yf cos- ži. )* 
yı ti-ra Ol 


EE r1,/801% -yf + ri/8o1 -Yi + Vi 801° 1" 801” -rQ* 
ri,/80 e-yi* ¥1,/ 80 2 r= 


/ 2 2 / 2 2 

Bor -yı +t 801 “1 
/ 2 2 / 2 2 
Soi “Yi 7 801 “Y1 


r,1n (69) 
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When equations (68) and (60) are evaluated for several ratios of 
gap to semispan, the results shown in figure 1(b) are obtained for the 
span loading on the wing-body combination. 


Pressure Distribution on Body 


The differentiation of APip (eq. (67)) yields 


OMPig 2105 sıtı i ] E(k) 
a (00 3090 -fa - FS [2 x a 


k sin Y, 
x: [917-01 [1 T Ao(¥,k2) 1)) | ze) - "ux (10) 


so that 


4) (x(a) - E(k) - 


ApN _ dm t(s*-ro*) (t#+r,7) 
z * g(g?-t")(s?t?-rg 


E ] at) (o t*K(ki) , x Xs mg) É " 


K(k) (ir E)? “o s( tr 2) 
As) | | E York) - | (72) 


If t is set equal to ro in equation (71), the loading coeffi- 
cient for the case of zero gap is obtained as 


EST o 


(72) 
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The loadings on a spanwise strip of the wing-body combination, 
obtained from equations (71) and (62), are shown graphically in figure 
2(b) for several ratios of gap to semispan. 


Lift of Body 


The lift of the body is found in the same manner as in the first 
problem, namely, by integrating the span loading over the diameter. The 
integration is carried out in the §, plane for convenience. When 
equation (68) is substituted into the equation 


Ti dip 
E o  Ayi x 


the expression for lift on the body is 


t3 = z (sox? =r?) | 1-2 aS | + ut L- | AdyY,k1) + 
B 


2ri?K(ki) K( kı) 
———— E - Oty 801-117 E(k) - = | 
Xt; /sq*-ri* sin 


rj*K(k er = i 
Adis) + E |. FEE sis) -E a 


tty 801 -r1 


with the moduli and arguments of the elliptic integrals evaluated for 


81 = 801» 


(73) 


The expression for the lift of the wing-body combination in the 
presence of a gap is obtained by adding equations (73) and (65) and is 
given by 


i -r —— — > A A m Am am A me al ys s T Ia gp A tpa n a m rt m e m 
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(5). = zb + (801% - r4?) | - 2 aay |} + Att) | > ari? + 


2t. /S012-r,? (200) - K(ki) + ui - 2 2 ks) (74) 


where 84 = So, is used throughout. 


As previously mentioned, no effect of the afterbody is included in 
the analysis of the lifting properties of the body, so equations (73) 
and (74) are strictly applicable only to wing-body combinations having 
no afterbody. 


It is interesting to note that the lift of the wing-body combination 
per unit angle of wing deflection is identical to the lift of the wing 
of the combination (exclusive of body nose) per unit angle of attack 
(cf. eqs. (74) and (35)). 


For the case of zero gap (ti = rı), equation (74) may be written 


Po 
redi ro Ne]? Bo 
Clara uS UE € | 1 - E T 


14 o fo 
2 
2 Bo -L Bo 
A Dc Yo? 
8o Bo 


where Cin = 5 AS = lift coefficient of a slender triangular wing. Equa- 
tion (75) agrees with a result obtained in reference 10. 
Figures 3(b) and h(b) present the variation of lift with gap width 


for the wing, body, and wing-body combination for ratios of body radius 
to wing semispan equal to 0.5 and 0.216. 


en —À— À ——— ee —— 
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EFFECTIVENESS OF ATI.-MOVABLE WING 


The effectiveness of a control surface in producing lift is often 
of interest to designers of aircraft. It is usually expressed in terms 
of an effectiveness parameter defined by 


acy, 

-2| 3 
6 dCr, 
da 


Ci] =constant 


For the wing-body combination investigated in this report, the 
variation of the effectiveness parameter with gap width is shown in 
figure 5. Figure 5(a) shows results obtained by omitting the lift of 
the nose, which according to theory is invariant with gap, and by omit- 
ting the lift of any afterbody, which cannot be calculated in the present 
analysis. In figure 5(b) the nose lift has been included in the calcula- . 
tions, so that the latter figure pertains to a slender wing-body com- 
bination having no afterbody. 


DISCUSSION OF RESULTS 


Application of slender-body theory to the determination of the Lift 
of a slender wing-body combination having & gap between wing and body 
indicates that comparatively large losses in lifting pressures, circula- 
tion, and lift are caused by even minute gaps; however, it is antici- 
pated that these large losses will not be realized in practice for very 
small gaps because of the viscous and possibly the compressible proper- 
ties of actual Pluids. Experimental investigation is necessary to deter- 
mine the range of gap dimensions and of angles of attack for which the 
simplified theory is valid. 


If attention is confined to the wing-body combination exclusive of 
nose and afterbody, it is evident from figures 3 and 4 that for a given 
gap width the percent loss in lift is greater when both wing and body 
are at an angle of attack than when only the wing is set at an angle of 
incidence. Comparison of figure 3 with figure 4 shows that the percent 
loss in lift due to gap increases with increasing ratios of body radius 
to wing semispan. Figure 5(a) indicates that the presence of a gap 
tends to reduce the difference between the lift of the wing-body com- 
bination due to angle of attack and that caused by wing deflection. 

The influence of the ratio of body radius to wing semispan is such that 
for larger ratios the effectiveness of the wing as a control surface is 
diminished. ' 
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The effect of the lift of the nose in the angle-of-attack case is 
to reduce the percent loss in lift of the complete configuration for a 
given gap width from that given for the portion of the configuration 
above; thus the effectiveness of the wing as a control surface is 
Lowered from the values given in figure 5(a) to those shown in figure 
5(b). It appears from figure 5(b) that, with no afterbody, the control 
effectiveness of the deflected wing may increase or decrease with 
increasing gap width accordingly as the ratio of body radius to wing 
semispan is small or large. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., May 12, 1954 
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APPENDIX A 


NOTATION 


8 
aspect ratio, l — 


lift coefficient, Y 
qS 


complete elliptic integral of the second kind with modulus k, 


1/2 
f /1 - k?gin^O do 
O 


incomplete elliptic integral of the second kind with argument 


y. 
Y and modulus x, f / 1 - k*sin*e de 
O 


incomplete elliptic integral of the first kind with argument 


do 


y 
Y and modulus 2 ———— mars mam 
O l - k?gin?o 


gap width between fuselage and wing panels, t - ro 


moduli of elliptic integrals, k = 


complementary modulus, /1 - X 
complete elliptic integral of the first kind with modulus k, 


[^ 25 
° 1 - ksinfo 


lift force 


slope of leading edge of wing panel, triangular wings, = 


free-stream Mach number 


ON 
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= loading coefficient, A 

q free-stream dynamic pressure, = Poo? 

fo radius of fuselage (cylindrical portion) 

8 local semispan of wing 

Bo maximum semispan of wing 

s*, gg* local and maximum semispan of wing, respectively, when gap is 
Zero 

5. combined areas of both wing panels 

t distance from body axis EN inside edge of wing panel 

u,V,W perturbation velocity components in the x,y,z directions, 
respectively 

Au, Av jump in velocities across the z = 0 plane (Uy - uj, Vu - Vi) 

Vr radial component of perturbation velocity in a plane perpen- 
dicular to the x axis 

Vo free-stream velocity 

X,y; Cartesian coordinates 

a angle of attack of body axis 

r total circulation about a wing section 

O wing deflection angle with respect to body axis 

€ apex angle of wing panel 

Z(W,Xx) Jacobi's zeta function, E(¥,k) - B F(v ,k) 

7 y” 

0 polar coordinate in a plane perpendicular to the x axis 

Ao(¥,k) Heuman's lambda function, 
2 fiso - K(k)] F(¥,k') + K(x)B( ys") } 

E complex variable, y + iz 


3h 


WB 
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polar distance in a plane perpendicular to the x axis 
free-stream density 
perturbation velocity potential 


jump in ỌPỌ across the z = O plane 


[aoe ae 
arguments of elliptic integrals, Y = sin”? a; 
B1 
27.2 2 2 
81 -t;* t? -y 
gin"i (s;*-yi4)ti*. Y, = sin"l 
(s e t =) ze 4 
1 “ti JJi 


Subscripts 


(81? -rf ) yf 


(817 -y3?)ri* 


o 
lI 


complex plane resulting from application of Joukowski trans- 
formation to the physical complex plane (except when used 
in the elliptic integrals) 


body 
lower surface of z = O plane or of wing panel 
upper surface of z = O plane or of wing panel 
wing 


wing-body combination 
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Figure 2.— Lifting pressure distribution on a sp 
a wing—body combination for se veral values of gap width. 
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Figure 3.- Variation with gap width of lift of a slender wing-body combination, g= 0.500 " 
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Figure 4.- Variation with gap width of lift of a slender wing-body combination, 
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Figure 5.- Variation of effectiveness of all-movable wing of a slender wing-body 
combination with gap size. 
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